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Abstract. The Gutzwiller trace formula provides a semiclassical approximation for 
the density of states of a quantum system in terms of classical periodic orbits. In 
its original form Gutzwiller derived the trace formula for quantum systems without 
spin. We will discuss the modifications that arise for quantum systems with both 
translational and spin degrees of freedom and which are either described by Pauli- or 
Dirac-Hamiltonians. In addition, spectral densities weighted by expectation values 
of observables will be considered. It turns out that in all cases the semiclassical 
approximation yields sums over periodic orbits of the translational motion. Spin 
contributes via weight factors that take a spin precession along the translational 
orbits into account. 



Thirty years ago, after several intermediate steps the Gutzwiller trace formula 
resulted from a detailed semiclassical investigation of the time evolution 
in quantum mechanics. It opened the way for the application of semiclassi- 
cal methods to many problems that were so far believed to lie beyond the 
capability of semiclassics. The most prominent example being semiclassical 
quantisation rules for classically non-integrable systems, for which Einstein 
[H already in 1917 had shown that the usual Bohr-Sommerfeld type quantisa- 
tion methods fail. Gutzwiller, however, devised a semiclassical expansion for 
the density of states in terms of a sum over the classical periodic orbits for a 
huge class of quantum systems, which in particular includes classically chaotic 
systems. Shortly afterwards, but seemingly independently, also mathemati- 
cians became interested in such trace formulae. They devised mathematical 
proofs for various versions of the trace formula, beginning with the work of 
Colin de Verdiere , and Duistermaat and Guillemin |^ . 

What was lacking so far, however, was a trace formula for quantum sys- 
tems with a priori non-classical degrees of freedom as, e.g., spin. In such cases 
it is not immediately clear what the corresponding classical system is whose 
periodic orbits enter the trace formula, and how the non-classical degrees of 
freedom have to be taken into account. Even for systems with a classically 
integrable translational part there do not exist Bohr-Sommerfeld (or EBK) 
type quantisation rules for the eigenvalues of a Dirac-Hamiltonian, although 
already in 1932 Pauli Q began to generalise the WKB method to the Dirac 
equation. Pauli's undertaking was only completed in 1963 by Rubinow and 
Keller [Q, and it took again some 30 years before Emmrich and Weinstein 
M proved that due to geometric obstructions for the Dirac equation EBK- 



2 Jens Bolte 



quantisation rules are generally impossible. It therefore seems to be especially 
desirable to have a Gutzwiller trace formula for Dirac-Hamiltonians available. 
Below we explain how one proceeds, if the Hamiltonian at hand describes a 
quantum system with a spin 1/2 coupled to the translational motion, be it 
a relativistic or a non-relativistic situation. The method that is used was de- 
veloped in Q , where most of the details can be found that cannot be given 
here. 



1 Dirac- and Pauli-Hamiltonians 



In the following we will consider relativistic and non-relativistic particles 
with mass m, charge e and spin 1/2 in external static electromagnetic fields 
E{x) = —W(p{x) and B{x) — V x A{x). In the relativistic case the quantum 
dynamics are generated by a Dirac-Hamiltonian 



A{x) + mc^fi + e(p{x) , 



(1) 



where a = (Q;i,Q;2,a3) and f3 = ao are hermitian 4x4 matrices satisfying 



the algebraic relations ttfj^a^ + OiuOLfi 
read 



2(5p,y. In Dirac representation they 



a- 
a- 
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ao = /? = 



l2 









(2) 



where cr = (ci, (72, ca) denotes the Pauli matrices and I2 is a 2 x 2 unit 
matrix. As it is well known, in leading non-relativistic order H^) is replaced 
by the Pauli-Hamiltonian 



1 
2m 



1 



Uix) 



■ e(p{x) 



2mc 



B{x) 



(3) 



For further information, see [Q. We remark that for the following the ex- 
plicit form of the right-most term in (H), describing a coupling of spin to the 
translational degrees of freedom, is inessential. It can be any operator of the 
form hC{j'V, x) ■ cr, where the components of C{p, x) are suitable functions 
on phase space; e.g., C{p, x) = ^^^^|^| jj!^^^ (x x p) would yield a spin-orbit 
coupling in a spherically symmetric potential. 

In applications one is sometimes also interested in describing a coupling 
of spin and translational degrees of freedom that is semiclassically strong. In 
this case one would consider a (Pauli-) Hamiltonian of the form 



Hp 



1 fh 

2m \ 1 



Uix) 



eip{x) 



D\ -V,x 



(4) 



where again the components Dk{p,x) are suitable functions on phase space 
(independent of h) . 
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1.1 Weyl Representation 



All of the above Hamiltonians can be represented as Weyl operators, i.e., their 
action on n-component spinors ip{x) = {ipi{x), . . . ,ifin{x))'^ , where n — 4 
applies to the case of Dirac spinors and n = 2 to Pauli spinors, can be given 
as 



1 



(27rft)3 



H (p, eiP<^-y^ 4,{y) dpdy 



(5) 



On the right-hand side the Weyl symbol -ff (p, x) is a function on phase space 
taking values in the hermitian n x n matrices. For the Hamiltonians (|l|),(H) 
and (Q) these arise upon replacing jV by p. 

In the case of the Dirac-Hamiltonian (|^) the symbol H£,{p, x) is a hermi- 
tian 4x4 matrix with the two doubly degenerate eigenvalues 



H (p, x) = eip{x) ± y/ {cp — eA{x)y + rrfid^ 



(6) 



These eigenvalues can be recognized as the classical relativistic Hamiltonians 
for spinless particles and anti-particles (— ). The fact that the correspond- 
ing eigenspaces are two dimensional reflects the quantum mechanical spin 1/2 
of the particles and anti-particles. More precisely, at each point (p, x) of phase 
space the eigenspaces of Hj:,{p, x) can be viewed as Hilbert spaces of a spin 
1/2, one for particles and one for anti-particles. See for further details. 

The situation is similar for the Pauli-Hamiltonian (^), except for the ab- 
sence of anti-particles: Hp{p,x) is a hermitian 2x2 matrix. 



Hp{p,x) = 



1 / e \ 2 
— [p- -A{x)) +eip{x) 
Zm \ c / 



I2 - 



eh 
2mc 



B{x) 



(7) 



whose classical part (independent of h) is proportional to I2. This principal 
symbol therefore has one doubly degenerate eigenvalue, Hq{p, x), which sim- 
ply is the factor multiplying I2. It can readily be identified as the classical 
non-relativistic Hamiltonian for spinless particles. As in the previous case 
spin is represented by the corresponding two dimensional eigenspace of the 
principal symbol, which in this case is trivial. 

The Hamiltonian (^), however, leads to a different interpretation of its 
symbol. 



Hp,{p,x) = 



J- (p- -Mx)] +etp{x) 
2m \ c ' 



I2 + -D(P, X) ■ (T 



(8) 



since this is independent of h and thus has to be considered as a classical 
quantity in its entirety. Its two eigenvalues 



H^IHp,x) = (p--A{x)\' + e^{x) ± \D{p,x) 
Zm V c / 



(9) 
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are classical non-rclativistic Hamiltonians of a particle with fixed spin up 
or down, respectively. Here 'up' and 'down' are defined with respect to the 
direction of D at the respective point (p, a;) in phase space. The eigenvalues 

are non-degenerate as long as they are different, i.e., away from mode 
conversion points where D{p, x) = 0. 

From the above discussions one may anticipate that in the semiclassical 
considerations of the Hamiltonians (|l|) and (^) to follow (see also |^) there 
will occur classical translational dynamics that are uninfluenced by spin, and 
quantum mechanical spin dynamics driven by the classical translational mo- 
tion. In contrast, in the case of the Hamiltonian the classical translational 
motion will depend on the (fixed) direction of the spin and there will be no 
additional spin dynamics. 

1.2 Spectra 

Gutzwiller's trace formula provides a semiclassical expansion for the quantum 
mechanical density of states. Since this quantity a priori requires a Hamil- 
tonian with a discrete spectrum, we now want to discuss the spectra of the 
Hamiltonians @), (|) and (|). 

Dirac-Hamiltonians typically possess continuous spectra. Indeed, if the 
electromagnetic fields vanish as \x\ — > oo, Hj^ has an essential spectrum 
consisting of two half axes, (— oo, — mc^] U [-I-toc-^, +oo), see [0 for precise 
statements. Thus, since we are interested in the discrete spectrum of a Hamil- 
tonian, we have to localise in energy to within the gap (— mc^, -|-mc^) of the 
essential spectrum. To this end one can choose a smooth function x(-E') that 
vanishes outside of some interval contained in (— mc^, +mc^) and then con- 
siders the Hamiltonian x{Hd)- This now has a purely discrete spectrum with 
eigenvalues x(-^n); if the EnS are the eigenvalues of Hd with |i?„| < rru? . 

The same procedure can be applied to the other Hamiltonians, if situa- 
tions arise where their spectra are not purely discrete or when one is only 
interested in certain spectral stretches. For the purpose of deriving a semi- 
classical trace formula one then considers a truncated time evolution operator 

U^{t)^e-TH^X{H) , (10) 
whose spectral expansion in position representation reads 

K^ix,y,t) = {x\U^{t)\y) =J2xiEn)M^)Myy e'^''''' , (11) 

n 

where il'nix) denotes the eigenspinor of H associated with £'„. 

The Gutzwiller trace formula has found many applications in the field of 
quantum chaos P,|ll[| in which the principal questions are associated with the 
distribution of eigenvalues and eigenfunctions of quantum Hamiltonians in 
relation to properties of the corresponding classical dynamics. In this context 
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one first confines oneself to some spectral interval / that contains Nj < oo 
eigenvalues; universal statistical properties then emerge in the semiclassical 
limit Ni — > oo. Here we implement this procedure by first choosing an interval 

I = [E - nuj,E + huj] , (12) 

and then performing the limit H ^ 0; and although the length of / shrinks to 
zero, Nj diverges in this limit (see below). In this approach the localisation 
in energy described above hence appears to be very natural. 



1.3 Observables 

The quantum mechanical observables we consider are (bounded) Weyl oper- 
ators, i.e., they can be represented as in (^). Their symbols B{p, x) are then 
functions on phase space taking values in the hermitian n x n matrices, and 
we suppose that they allow for asymptotic expansions in h, 

B{p,x) ^Y.^'' Bk{p,x) . (13) 

fc>0 

The precise meaning of such expansions is explained in fl^ . The ft-indepen- 
dent term Bo{p,x), the principal symbol of the operator B, represents the 
classical observable associated with i?, at least concerning the translational 
degrees of freedom. Expectation values of observables can be represented in 
terms of the symbol once one introduces a matrix valued Wigner transform 
of a spinor ip{x) — (tpi^x), . . . ,'0„(a;))^ through 

W[^p]kl{p, x) := J c-T^P-y~k{x - \v) + ^y) dy . (14) 
Then 



{4',BiP) = j^-^ J J tr (Wmp,x)B{p,x)) dpdx . (15) 
2 Trace Formula 

Gutzwiller's approach to the trace formula was to depart from a semi- 
classical expansion of the time evolution operator in position representation 
(pr]). Expressing this in terms of a Feynman path integral and evaluating it 
in leading semiclassical order with the method of stationary phase, he arrived 
at a representation of the kernel K{x,y,t) in terms of a sum over the clas- 
sical trajectories connecting y and x in time t. In this context he made the 
important observation [ p^ that for not too small times t each term in this 
sum must contain an extra phase factor that essentially consists of the Morse 
index of this trajectory. 
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Here we are interested in a trace formula for quantum systems with spin 
that yields a semiclassical expansion of a weighted and smeared spectral 
density, where the weights are provided by the expectation values of an ob- 
servable B in the eigenstates of the Hamiltonian and the smearing ensures 
convergence of the sums involved, see also [Q. To this end the 'spectral' side 
of the trace formula is given by 

J2 X{En) ii^n, Hn) P = Tr- / ~p{t) 6*^' B U^{t) dt . 

In this expression p{E) is a smooth, regularising function with Fourier trans- 
form p(t) that is also smooth and vanishes outside of some compact interval. 
The trace on the right-hand side will now be calculated in position represen- 
tation, after the leading semiclassical order of ([III), i.e., the appropriate Van 
Vleck-Gutzwiller propagator, has been introduced. 



2.1 Van Vleck-Gutzwiller Propagator 

In an alternative approach to Gutzwiller's semiclassical expansion of the prop- 
agator ( pi] ) , making use of a Feynman path integral, one represents the ker- 
nel in a way that is closely related to the WKB method. More precisely, one 
chooses the ansatz 

K{x,y,t) = -i— /E f-V et(^^('>'-S.*)"«€) d^ , 

^^'"^> J k>o^'^ (17) 

and determines the phase S{x,^,t) and the coefficients ak{x,y,t,^) of the 
amplitude by requiring that K{x,y,t) solves the Schrodinger (Dirac, Pauli) 
equation to arbitrary powers in h with an appropriate initial condition at t = 
0. In the cases of the quantum Hamiltonians (|l]), (^) and (jj) it turns out that 
the corresponding phases S{x,^,t) have to be solutions of Hamilton- Jacobi 
equations with H^(n,x), Hq{p,x) and H'^^^{p,x), respectively, as classical 
Hamiltonians (see (^-(|^)). The matrix valued coefficients ak{x, y, t, ^) of the 
amplitudes are determined by a hierarchy of (transport) equations that can be 
solved order by order in h, starting with the leading expression ao(a;, y, t, ^), 
which apart from the contribution of the translational motion contains the 
leading order of the spin dynamics. Knowing the solutions S{x,^,t) and 
ao(a;, y, i, ^), one then still has to evaluate the integral ( p7| ) over ^ with the 
method of stationary phase. The details of this calculation can be found in 

Since S{x,^,t) is a solution of a Hamilton- Jacobi equation, it is a gen- 
erating function of a canonical transformation, {p,x) i-^ (^i^;). Here p = 
\/xS{x, ^,t) and z — V^S{x,^,t), such that {^,z) and {p,x) are starting 
and end points, respectively, of a solution of the equations of motion gener- 
ated by the respective classical Hamiltonian. The stationary points of the 
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phase in (Jl^) are now uniquely related to classical trajectories from y to a;, 
since the condition of stationarity reads y = V^S{x,^^^^,t). The stationary 
point itself hence is the momentum of the associated trajectory at time 
zero. Thus the sum over the stationary points leads to a sum over classical 
trajectories; see |^ for the explicit form of these sums. In the cases of the 
Dirac- and Pauli-Hamiltonians (|l|) and (^) the trajectories that contribute 
are purely translational without contributions from spin (since the respective 
classical Hamiltonians contain no spin). The latter only contributes through 
weight factors that come from the leading term ao(a;, y, ^g^) of the ampli- 
tude. Since this is evaluated at the stationary point ^g^, the spin contribution 
derives from a precession along the associated classical trajectory. For a fur- 
ther discussion see also the next subsection; more details are described in 

Only in the case of the Hamiltonian (^) is the spin already directly con- 
tained on a classical level, see (^). Here one has to deal with the (two) trans- 
lational dynamics of a particle whose spin is tight to the direction of the 
external 'field' D, such that it follows this direction adiabatically along the 
trajectories of the particle. This adiabatic motion of the spin is reflected in 
the occurrence of a certain geometric phase, which has in the context of a 
Bohr-Sommerfeld quantisation been introduced by Littlejohn and Flynn JTsf , 
see also [IGlOl for a discussion in the context of a trace formula. 



2.2 Semiclassical Spin Transport 

In the cases of the Dirac- and Pauli-Hamiltonians (|l|) and (||) there indeed 
is a dynamics of the spin degrees of freedom, which is driven by the classi- 
cal translational motion. These driven dynamics derive from the transport 
equation for the lowest order amplitude ao{x,y,t,^^^) after separation of 
the purely translational part. The spin transport equation that hence results 
reads [|j 

d{p, x,t)+i C(p{t),x{t)) ■ (T d{p, a;, i) = , (18) 

with initial condition d{p,x,0) = I2. Its solution, the spin transport ma- 
trix d{p,x,t) e SU(2), propagates the (quantum) spin 1/2 along the clas- 
sical trajectory {p{t),x{t)) starting at {p,x). The vector C{p,x) depends 
on which Hamiltonian one considers and contains the fields E and B; for a 
Pauli-Hamiltonian Hp the quantity C is precisely the one described below 
(||). Moreover, if the vector s{t) denotes the expectation value of the (nor- 
malised) spin operator cr in a two-component spinor u{t) = d{p,x,t) u{0), 
this 'classical' spin obeys the equation 

sit)^C{p{t),x{t))xsit) (19) 

of classical spin precession and thus provides an Ehrenfest relation for the 
spin. In the relativistic case (foh yields the well known Thomas precession, 
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which Rubinow and KeUer j^] were the first to derive semiclassically from the 
Dirac equation. 

One can now combine the Hamihonian translational motion and the spin 
dynamics that are driven by the former one into a single dynamical system 
on a combined phase space, (p(0), a?(0), s(0)) ^ {p{t),x{t),s{t)). In ergodic 
theory such combinations are known as 'skew products' of the two types of 
dynamics. In applications to quantum chaos the ergodic properties of pre- 
cisely these combined dynamics determine the 'quantum chaotic' properties 
of the quantum system, see |l7|,^ 19 1. 



2.3 Semiclassical Trace Formula 

As mentioned earlier, the 'semiclassical side' of the trace formula emerges 
upon introducing the leading semiclassical order of the propagator ( [l7| ) on 
the right-hand side of ( [l6| ) and evaluating all integrals, which involve the 
variables {^,x,t), with the method of stationary phase. Since this method 
requires all stationary points (^g^, a^st, ^st) to be non-degenerate, the trace 
formula can only be derived under appropriate conditions on the classical sys- 
tems. In particular, since the stationary points are such that the phase space 
points (^st, Xst) lie on a periodic orbit with energy E (that appears on the left- 
hand side of (|l6|)) and period tst, the conditions have indeed to be imposed 
on the periodic orbits. One such condition is that E must not be a critical 
value of the relevant classical Hamiltonians, i.e., {VpH{p, x), VxH{p, x)) ^ 
for all {p,x) on the energy shell He = {{p,x); H{p,x) = E}. This condi- 
tion ensures that all stationary points with tgt = are non-degenerate. The 
corresponding points (^gj,£Cst) make up all of the energy shell He and the 
contribution of these stationary points yields the leading semiclassical term 
(also called Weyl term) on the right-hand side of the trace formula. 

For the following we restrict our attention to the case where all non-trivial 
periodic orbits 7 (i.e., T-y = tst 7^ 0) are isolated and non-degenerate. This 
means that their monodromy matrices Mj, describing the linear stability of 
the orbits, have no eigenvalues one. This does not exclude elliptic orbits, if 
these are isolated, but only parabolic (marginally stable) ones. Furthermore, 
we give the trace formula explicitly for the case of a Dirac-Hamiltonian (Q) 
and with the inclusion of an observable B. The other cases follow from this 
trace formula by specialising to the appropriate simplified situations; e.g., 
the Pauli-Hamiltonian (^ has no contribution from anti-particles and there- 
fore only contains one type of (translational) classical dynamics. In contrast, 
the Hamiltonian (^) does lead to two types of classical dynamics, generated 
by (^), but has no independent spin dynamics; there only is an additional 
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geometric phase as described in [ p^jjl^ j^ . Now the trace formula reads |p|,p^ 

E X{En) (^n, Hn) P (^^^) (20) 

P(0) xiE) ( 2 ^^-^E ^ ^ ^^^^^ trBo"^) + 0{h-^) 



+ E ^(^) ^ trs;;^" e*^^- (1 + 0(?.)) . 

The first term on the right-hand side is the so-called Weyl term and yields 
the leading semiclassical approximation of the left-hand side. It essentially 

—tE 

contains the averages of the principal symbol Bo{p, x) of the observable 
over the two energy shells f2^ . The two contributions are weighted according 
to the relative volumes of the respective energy shells. The characteristic 
quantities appearing in the sum over the two types 7^ of periodic orbits are 
the same as in the case of the Gutzwiller trace formula for the spectral density 
of a quantum system without spin: there is an exponential factor with the 
actions S-y{E), and an amplitude that contains the well known part 

= |detV7-l)l^/' ' ^^^^ 

with the associated primitive period, the Maslov index /i^, and the stability 

denominator. The additional factors come from the regularisation and from 

^± 

the presence of the observable. In particular, Bq denotes an average of 
the projection of Bo{p,x) to one of the eigenspaces corresponding to the 
eigenvalues (^, weighted with the spin transport matrix d{p, x, t), along the 
periodic orbit 7^. Details of this average are described in jl^ . 

A trace formula for a truncated spectral density, with only contributions 
from the discrete spectrum of Hd, can be obtained from (^ by choosing the 
observable i? = id and removing the regularisation provided by the function 
p. The result is M 



dx(i?) = E^(^")'5(^«"^) 



-^(^^ '^"'i.'.r''' ^^^^'^^ (22) 

+ 27r?i ^ |det(M^±-l)|i/2^ U + ^WJ , 

where here d^± denotes the spin transport matrix associated with the trans- 
port of spin once along the periodic orbit. One immediately observes that in 
contrast to the Gutzwiller trace formula without spin two types of classical 
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dynamics (particles and anti-particles) occur, and spin contributes a factor 
of two in the Weyl term as well as weights trd^i for the periodic orbits 
describing the effect of spin transport. 

Obviously, the Weyl term in ( p2| ) yields the semiclassically leading contri- 
bution to the spectral density. This allows to derive the leading semiclassical 
asymptotics of the number iV/ of eigenvalues in the interval (|l2|), which is 
given by 2hLu times the Weyl term, and hence Nj diverges as h^^ in the 
semiclassical limit. 



3 Applications 

Many applications of the Gutzwiller trace formula originate from problems 
in the field of quantum chaos (see, e.g., p|jllf| ). Apart from the question for 
semiclassical quantisation rules (see, e.g., ^0|) one of the major successes in 
this field was Berry's semiclassical analysis of spectral two-point correlations 
[ pT| based on the trace formula. In a certain range of validity, which stems 
from the so-called diagonal approximation Berry employed, he verified that 
the two-point correlations of energy levels of classically chaotic quantum sys- 
tems (without spin) follow the predictions of random matrix theory (RMT). 
Subsequently, Bogomolny and Keating |Q extended Berry's result in that 
they went one step beyond the diagonal approximation. 

The first application of trace formula techniques, without, however, having 
a complete trace formula available, to quantum systems with spin goes back 
to Frisk and Guhr . They considered a quantum Hamiltonian of the type 
(^ describing spin-orbit coupling in certain billiards. In spirit, they applied 
the trace formula in 'reverse direction' in that they used quantum energy 
levels in order to obtain information about the contribution of various types 
of periodic orbits. Similar studies can be found in 

A second type of applications concerns an extension of the semiclassical 
analysis of spectral two-point correlations to quantum systems with spin 1/2. 
If a time-reversal symmetry is present, Kramers' degeneracy implies that all 
energy levels of a quantum system with half-integer spin are doubly degen- 
erate. After removal of this systematic multiplicity, the spectral statistics 
should be described by the Gaussian symplectic ensemble (GSE) of RMT, if 
the corresponding classical system is chaotic. Without time-reversal symme- 
try Kramers' degeneracy is absent, and the relevant ensemble of RMT is the 
Gaussian unitary one (GUE). In [0 Berry's semiclassical approach is carried 
out with the trace formula ([20|), or (p2[), as a basis. It is shown that, within 
the same range of validity as in pl| , indeed the two-point correlations agree 
with the GSE or GUE, respectively. Moreover, in the Bogomolny-Keating 
method is carried over to the case of spin 1/2, with the same findings. In both 
studies, apart from a chaotic translational motion, ergodic properties of the 
combined dynamics described below (|9|) are needed in order to obtain an 
agreement with the RMT predictions. 
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The Weyl term of the trace formula (|2^) moreover allows to determine 
a semiclassical average of the expectation values {ipn, Btpn) in eigenstates of 
Hd with E„ e I. The result is given by Q 

1 , , - , , 1 vol/2+trS^^ + vol/2p trB^^ , , 

™ aT E ^^") = 9 1° + ^ ^o~ ' ^2^) 

?i-^0 iV/ 2 VOliiX + VOlJ/p 

E„£l ^ ^ 

and does not require the classical system to be chaotic. The only require- 
ment is that the set of periodic orbits on both energy shells must be of 
measure zero, which is a comparatively weak condition. Only if one requires 
individual expectation values to approach the expression on the right-hand 



side of (23) as ?i ^ one needs a stronger condition on the classical side, 
since then no cancelations on the left-hand side are allowed. For the case of 
Pauli-Haniiltonians (||) it has been proven ||l8| (see also |l9|) that almost all 
expectation values indeed converge to the equivalent of the right-hand side of 
(p3|), if the combined dynamics of translational and spin degrees of freedom 
are ergodic. For quantum systems without spin such a result had been known 
before under the notion of quantum ergodicity (see, e.g., ]25[|). 
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